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ABSTRACT 


In order to expand the predictive capability of single-point turbulence closure models to account 
for the early-stage transition regime, a methodology for the formulation and calibration of model 
equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is 
presented. First the decay of laminar disturbances and turbulence in mean shear-free flows is 
studied. In laminar flows, such disturbances are linear superpositions of modes governed by 
the Orr-Sommerfeld equation. In turbulent flows, disturbances are described through transport 
equations for representative mean quantities. The link between a description based on a de- 
terministic evolution equation and a probability based mean transport equation is established. 
Because an uncertainty in initial conditions exists in the laminar as well as the turbulent regime, 
a probability distribution must be defined even in the laminar case. Using this probability dis- 
tribution, it is shown that the exponential decay of the linear modes in the laminar regime can 
be related to a power law decay of both the (ensemble) mean disturbance kinetic energy and 
the dissipation rate. The evolution of these mean disturbance quantities is then described by 
transport equations similar to those for the corresponding turbulent decaying flow. Second, ho- 
mogeneous shear flow, where disturbances can be described by rapid distortion theory (RDT), 
is studied. The relationship between RDT and linear stability theory is exploited in order to 
obtain a closed set of modeled equations. The linear disturbance equations are solved directly 
so that the numerical simulation yields a database from which the closure coefficients in the 
ensemble-averaged disturbance equations can be determined. 
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1 Introduction 


Demands on the range of applicability of turbulence modeling are increasing, and with these 
increasing demands has come the need to develop transition models within the context of the 
traditional Reynolds-averaged turbulence modeling. Such transition models would encompass 
the laminar/ transitional regimes, along with the turbulent regime, describing averaged flow 
properties such as the mean disturbance energy and dissipation rate. In order to develop a 
transition/turbulence model for wall-bounded flows, which uses the strategy of employing an 
intermittency function that interpolates between the laminar regime with its linear disturbances 
and the fully turbulent regime with its stochastic fluctuations, a consistent mathematical de- 
scription of the two regimes is necessary. 

In all flow regimes, disturbances are defined here as deviations from the ensemble mean. The 
laminar regime is defined as the region of the flow in which the ensemble mean velocity (zero in 
the case of the homogeneous flows considered in this study) corresponds to a stationary solution 
of the Navier-Stokes equation, and disturbances from this mean velocity are small enough in 
amplitude that their nonlinear interactions can be neglected. In this regime, the evolution of the 
disturbances is completely predictable from their initial state. Traditionally in laminar stability 
theory, disturbance fields are studied through the linear Orr-Sommerfeld equation, which de- 
scribes the evolution of individual infinitesimal disturbance modes. Even when a quantity, such 
as the disturbance energy, is studied, it is the evolution of the instantaneous quantity rather 
than an ensemble average that is investigated. In contrast, the turbulent regime is defined as the 
region where the flow is subject to stochastic fluctuations, arising from nonlinear interactions, 
which render the behavior of the disturbances unpredictable. In this case the disturbance field is 
studied through (modeled) transport equations which describe the evolution of mean turbulent 
correlations. The purpose of this study is to reconcile these two apparently disparate approaches 
through a common mathematical analysis. 

While the Orr-Sommerfeld equation accounts for viscous effects on disturbances, it neglects 
the influence of any nonlinear interactions. The behavior of finite-amplitude disturbances can 
be different and is more complex since nonlinear interactions occur. However, the intent here 
is to analyze, within the framework of (ensemble) mean disturbance transport equations, the 
behavior of the linear disturbance fields which characterize the beginning of a transition pro- 
cess, and ultimately lead to a fully turbulent field. The approach taken here is based on the 
observation that even in the laminar regime every flow is subject to an inevitable uncertainty 
in initial conditions. Therefore, although each individual disturbance evolves deterministically, 
a probability distribution must be introduced for the calculation of ensemble mean properties. 
This approach is similar to rapid distortion theory (RDT) in that it is based on linearized dis- 
turbance equations; however, the physical interpretation is different. RDT considers flows in 
which turbulence is fully developed and uses linearized equations to study the behavior of the 
disturbances under rapid distortion. RDT is usually applied to short time evolution and the 
effect of viscosity is neglected. In the approach taken here, there is no limitation on the period 
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of time evolution as long as the disturbances remain small, and the effect of viscosity is essential. 
Thus linear disturbances are considered in the early stages of transition where viscous effects 
must be taken into account. 

More recent extensions of RDT by Salhi, Cambon and Speziale (1997) have also exploited the 
connection with linear stability theory. They studied quadratic flows in a rotating frame to gain 
better insight into the dynamics so that a generalized stability criterion applicable to turbulent 
flows could be developed. They also considered the effect on single-point closure modeling - 
specifically the deficiencies in predicting elliptic flows. While the mathematical framework is 
similar in this study, the region of interest here is the early-stage transition regime. Nevertheless, 
this commonality further substantiates the basic assumption that a mathematical framework can 
be developed which will provide a set of transport equations capable of describing the flow (in 
a mean statistical sense) in the early-stage transition regime. 

First the decay of laminar disturbances and turbulence in mean shear-free flows is studied. 
The study of homogeneous turbulence has become an essential element in the calibration of 
turbulence closure models. The simplest of the homogeneous flows is decaying, isotropic turbu- 
lence, which traditionally has been used to determine the destruction coefficient in the modeled 
dissipation rate equation used in two-equation and higher-order closure models. 

In the second part the case of mean homogeneous shear is studied. Homogeneous shear flow is 
commonly used as a calibration flow for turbulence models because both turbulent transport and 
viscous terms can be removed from the transport equations for the turbulent correlations. An 
extensive amount of experimental (Tavoularis and Corrsin, 1981, Tavoularis and Karnik,1989) 
and numerical simulation data (Rogers and Moin, 1987) at low and moderate shear also exists 
which further aids in its role as a turbulent calibration flow. 

The purpose of this entire study is to use the solution of the disturbance evolution equations 
for mean shear-free and mean homogeneous shear flows as a database in the calibration of the 
evolution equations for the ensemble-averaged disturbance correlations. As a beginning, the 
focus is on a simple disturbance kinetic energy and disturbance dissipation rate (two-equation) 
description. In such a two-equation description of homogeneous shear, the terms in the kinetic 
energy equation are exact and require no modeling; whereas, in the dissipation rate equation 
both the production-of-dissipation and destruction-of-dissipation terms require modeling. The 
two closure coefficients associated with these terms are determined from the analysis presented 
here. Utilization of the disturbance evolution results as a reliable database is supported by the 
DNS results of Lee, Kim, and Moin (1990) who studied homogeneous shear flow at a high-shear 
rate. They showed that RDT results compared very well with the simulation results over the time 
period examined. The results herein also show good agreement with the DNS results (Lee, et al } 
1990), and are found to apply at much later times due to the energetic decay of the disturbance 
field in the parameter range studied. Thus, this database will be used to provide insight into 
the asymptotic behavior of important dynamic variables, as well as to provide the necessary 
information for the closure model calibration. The resulting closed disturbance dissipation rate 
equation can then be used in the formulation of a transition-sensitized turbulence model. 
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2 Boundary-Free Disturbance Fields 


Consider first the case of disturbances in a laminar, zero mean-shear flow with no boundaries. 
Solid boundaries do indeed play a key role in the dynamics of any laminar disturbances and 
the transition process itself; however, in order to establish the mathematical framework, this 
first example will be unbounded, as is the decaying homogeneous turbulent flow. For this 
initial-value problem, an arbitrary solution of the Orr-Sommerfeld equation (here, the linearized 
Navier- Stokes equations) is given by 


u 


i(x,i) = J d 3 ku°(k)e^ k x “ uk2t \ 


( 1 ) 


where k = (ki , it 2 , k%) is the wavenumber vector in the coordinate directions x = ( x,y,z ), 
respectively, and u is the kinematic viscosity. Because the initial conditions are uncertain, an 
ensemble of such disturbances is considered for which a probability distribution can be ascribed 
to the initial mode amplitudes, u°(k). The mean of this distribution is zero and the covariance 
is ^u°(k)u°(k')). The corresponding two-point correlation function of the disturbance field is 

then ,2 

(Ui(x, f)tij(x',»)) = / d 3 W 3k-( u? (k) U »( k ')) e '~>e-^ +l X 

where x' = x + r, and the homogeneity of {rq(x, £)tij(x',<)) implies that 

(«?( k)u° j (k')) = 5 3 (k + k')f ij (k). 

The two-point correlation function is then given in the form 


( 2 ) 


(3) 


Rij(r,t) = I d 3 k/ij(k)e' kr ‘ 


-2i/ k 2 t 


(4) 


with the corresponding energy spectrum tensor 

Eij(k,t) = f tJ (k)e~ 2l/k2t , (5) 

which is assumed to be analytic at the origin (Hinze, 1975). This yields a wavenumber distri- 
bution for f x j which satisfies both isotropy and incompressibility, 

f ij (k) = (k 2 6 i j-k l k ] )f(k 2 ), (6) 

where f{k 2 ) is a nonsingular, scalar function which can be related to the mean disturbance 
kinetic energy through 

K(t ) = \J d3 k Ea{k,t) = J d 3 k f{k 2 )k 2 e- 2vkH . (7) 

In general, f{k 2 ) must go to zero sufficiently rapidly at infinity so that the integral in (7) is finite 
at t = 0. A natural choice for f{k 2 ) is the normal distribution exp(-afc 2 ). With this choice of 
f(k 2 ), Eq. (7) then gives the mean disturbance kinetic energy 

K{t) = K 0 ^1 + 1 , ( 8 ) 
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where Kq is the initial value of K . Recall that for decaying homogeneous turbulence, the power 
law behavior for the final period of decay (Hinze, 1975; Batchelor and Proudman, 1956), is 
K{t) oc t" 5 / 2 . Thus, the decaying linear disturbances have the same temporal character as 
turbulence in the final period of decay. 

The covariance of the initial mode amplitudes given in (3) can then be written explicitly as 
(ti?(k)tt}(k')) = § 3 {k + k') (Sij - K 0 V( k), (9) 


v{k) = I (i) (10) 

is a probability distribution for the disturbance second-moments, normalized so that / d 3 kP(k) = 
1, and a is related to the variance of the distribution. With this distribution 'P(k), other second- 
moments can be calculated. For example, the ensemble mean energy dissipation rate 


/ dui f chii_ duj \ \ 

\dxj \dxj dx t J I 


( 11 ) 


can be written as 


e 


-v J d 3 kd 3 kV (k + k ')- x e -l,(fc2+ k ' 2)t [(uj(k)u-’(k')) fcjfcj + («?(k)u°(k')} k 3 k\ 
2uKq J d 3 k k 2 e~ 2ukH V{k) = ^-K 0 ^1 + ^) 


( 12 ) 


which, for this decaying flow, is simply given by the time derivative of the disturbance kinetic 
energy in (8), that is, K = —e. 

It is well known that for any power law decay of the mean disturbance kinetic energy, with 
power -p, the mean disturbance dissipation rate equation 


£ = ~Ce2J£. ( 13 ) 

is satisfied exactly, and the coefficient C £ 2 is given by 

C t 2 = 1 + (14) 

P 


For this case of mean, laminar shear-free flow with no boundaries, the coefficient C e 2 is then 1.40. 
In the turbulent case, values in the range of 1.80 — 2.00 have been deduced from experiments 
(Batchelor and Townsend, 1948; Comte-Bellot and Corrsin, 1966; Comte-Bellot and Corrsin, 
1971), although the values more commonly used are delimited by 1.83 — 1.92. 

In the turbulence case, there has been a considerable amount of research associated with the 
proper choice of decay rate (Frisch, 1995; Lesieur, 1995). The values obtained have ranged from 
a power law decay with exponent -6/5 (corresponding to a k 2 low wavenumber behavior of the 
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energy spectrum E(k) oc k 2 Ea( k)) to a power law decay with exponent -10/7 (Kolmogorov 
decay law corresponding to a A; 4 low wavenumber behavior). In the laminar disturbance case, 
the temporal and the wavenumber dependence of the energy spectrum tensor are both given 
explicitly from the disturbance mode solutions of the Orr-Sommerfeld equation together with 
a probability distribution. The probability distribution enforces analyticity at k = 0 and goes 
to zero sufficiently rapidly at infinity so that the integrals exist. A simple integration over 
wavenumbers then yields the decay law of the disturbance kinetic energy with exponent -5/2, 
in agreement with the final period of decay for the fully turbulent case. This agreement is not 
surprising, because in the final period of decay, viscous effects are dominant, with the small scale 
(high wavenumber) turbulence already decayed away. 

3 Wall-Bounded Disturbance Fields 


While it is encouraging to have identified a relationship between the decaying linear distur- 
bances and decaying turbulence in the boundary-free flow, a wall-bounded flow is more relevant 
for an analysis of linear disturbances which evolve through a transition stage and then into 
full turbulence. Shear-free decaying turbulent flows with boundaries have been the focus of 
experimental, (Uzkan and Reynolds, 1967; Thomas and Hancock, 1977) theoretical, (Hunt and 
Graham, 1978) and numerical simulation (Perot and Moin, 1995) studies. The motivation for 
these studies was to assess the inhibiting effect on the turbulent fluctuations due to the presence 
of the wall in the absence of any mean shear arising from the relative motion between the mean 
flow and the wall. Here, the corresponding case of disturbances in a laminar, zero mean-shear 
flow with a boundary will be studied. For laminar disturbances, the wall has the effect of fixing 
the phase of the Orr-Sommerfeld modes in the wall-normal direction. 

The wall- normal direction is the y direction, and the wall boundary is located in the {x,z) 
plane at y = 0. In this flow, the modes are given by 

(15) 

where now k = (fci,0, k$) is the wavenumber vector in the coordinate directions (x, y,z), respec- 
tively, and oj is the complex frequency of the disturbances. An analysis (Grosch and Salwen, 
1978) of the Orr-Sommerfeld equation shows that bounded solutions only exist when to is pure 
imaginary (u? = i = y/—l) and 

u t = —(1 + A )vk 2 < — vk 2 , (16) 

where the parameter A(> 0) has been introduced. For three-dimensional disturbances, Squire’s 
transformation is used, and the Orr-Sommerfeld equation determines the velocity disturbances 
u 2 and k\u\ + k$u 3. The individual disturbance components, *1 and U3, can then be determined 
in one of two ways. First, the pressure field can be determined from the linearized 112 momentum 
equation and then substituted into both the u\ and U3 linearized momentum equations. Second, 
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the kinematic equation for the y component of vorticity i(k^u\ — k\u 3) can be coupled with the 
relation between 112 and (k\ui + ^3^3) from the Orr-Sommerfeld equation to determine u\ and U3 
individually. Because the solution for ( k\U\ + k^u^) and 112 obtained from the Orr-Sommerfeld 
equation is an irrotational wave of arbitrary amplitude, and the solution of the vorticity equation 
has an independent amplitude, it is necessary to introduce a new parameter, /r, for the ratio of 
these two amplitudes. The component disturbance modes are then given by 


tii(y,k) = i 

y cos i'/Xky) + (v^Ay - My) sin (v^Xfct/) - ye ky 

, (17) 

U2{y,k) - 

- cos (VX ky) H — y= sin (V\ky) + e ky 
v A J 

1 

(18) 

Ui(y,k) = i 

y cos (VJ-ky) + (v^Ay + My) s i n (>/Afcy) ye ky 

, (19) 


where k — \J k\ + k^ because fc-2 = 0. At the wall where the no-slip condition is applied, U{ 0, 
and far from the wall in the freestream, where e~ ky -> 0, the modes ii x are periodic functions 
of y. It should be recognized that these expressions for the velocity components, because of the 
no-slip wall boundary conditions, have fixed the phase of the disturbances in the y-direction. In 
the absence of the wall boundary condition, the exponential terms would not appear and the 

phase would remain arbitrary in this direction. 

The component disturbance modes (17) - (19) form a complete set (Salwen and Grosch, 1981) 
of non- normalizable eigenfunctions. A disturbance in the laminar regime is a linear combination 
of these eigenfunctions 

Ui(x,t) = fdX dn i 2 k <&(A, /.,*!, **)%(!>, k)el* kx ' “ + (20) 

which is normalizable in accordance with the finiteness of the disturbance kinetic energy. In 
Eq. (20), the initial mode amplitude <&(A, /x, &i, £3) is again an element of an ensemble with 
mean zero and covariance {^(A, (i, Aq, ^3) <J>( A', //, k\ k')). Since the modes u t are complex, the 
two-point correlation function is the real part of 

(n*(x, t)uj(x\ £)) = Jd\ dy d\' dy d 2 \<i d 2 k f {$ (A,/x, Aq, ks)<f>(\ ,y ,/c , A; )) 

■iffy, k)jij(y', k')e“ k ' *' - k "‘>e-'« 1 + + <‘ + (21) 

where * quantities denote complex conjugates. Homogeneity in the X\ and X3 directions then 
implies that ($*(A , /1, *1, *s)®(A', fc', * r )> contains 6 (k x - k') 5 (k z - k')V{k u k 3 ) where, due to 
invariance under rotations about the y-axis, V{k\, ks) is a function of A: 2 . It is still necessary to fix 
the dependence of the covariance (<h(A fcx , ^3)^** (A f , /i 7 , k',k')) on (A, /z, A',/x')- This is done by 
imposing an appropriate isotropy condition in the freestream, where the decaying exponentials 
can be set equal to zero in the component disturbance mode equations (17) - (19). In the 
boundary-free case, the disturbance field is isotropic, and the ensemble mean square of the three 
velocity components and of the three vorticity components are equal at each point throughout 
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the domain. In the bounded case, the wall fixes the phases of the laminar disturbances in the 
y direction even for large values of y, and it is not possible to impose pointwise isotropy as 
a condition on an ensemble of such disturbances. However, for large values of y (freestream), 
a weaker form of isotropy can be imposed on the disturbance field by averaging over the y 
direction. Specifically, this requires that in the freestream the ensemble mean square of the 
three velocity components and of the three vorticity components, integrated over y, be equal. 
The simplifying assumption must also be made that ($(A, y, k\, k 3 )$*(X', y' , k')) is diagonal in 
A and y so that 

<$(A, y, kuk^Q'iX', y\ k')) = <5(A' - \)5{y' - y)S 2 { k - k!)V{k x , k 3 ) p{ A, y) (22) 


where the function p(X,y) remains to be determined from the isotropy conditions. Taking 
yi sufficiently large that e~ ky « 0 for y > y x , one obtains from (17) - (19) 


lim — 
L—tco L 


dy <J|ui(x,t)| 2 J> 


1 fV l+L , 

lim y / dy (\u 2 {x,t)\ 2 ) 

L > — ^OO ]_j J yi ' 

lim y f V dyl\u 3 {x,t)\ 2 ) 

L — ► oc J y\ ' 


J dk\dk 3 dXdyV{k \ , k 3 )p(X, y) 


1 , \ \ 2 ^3 
2 (1 + A) fc2 +/i 


,,-2(1 +\)vk 2 t 


J dk\dk$d\diiP{k\, k$)p(\, y) 
J dk\dkzd\dfiP{k \ , k^)p(X 3 y) 


a + A)| + M 2 B 


:> -2(l + A )i/k 2 t 


(23) 


(24) 


(25) 


where the integration over y has eliminated the cross terms cos(\/AA;y) sin(\/A ky), while sin 2 (\/A ky) 
and cos 2 (v^A ky) have each yielded an average of 1/2. The terms containing k\k 3 have also van- 
ished since V{k x ,k 3 ) is even in k x and in k 3 . According to the assumed isotropy, the function 
p(X,y) must have the property that the right hand sides of Equations (23)-(25) must be equal. 
The condition that averaged over y, (|ui| 2 ) = (|^3| 2 ) is satisfied automatically, while the con- 
dition (|ui| 2 ) = (|u 2 | 2 ) yields the relation 

i(l + A + M 2 ) = l + {, ( 26 ) 

where use has been made of the identity 

J d 2 kk 2 F(k 2 ) = j d 2 kk 2 F(k 2 ) = l -J d 2 kk 2 F(k 2 ) (27) 

when F is any function of k 2 for which the integral exists. Calculating the vorticity Uj(x,t) from 
the modes u,(y, k)e (ik x - (i+^M 2 ^ with e~ ky = 0, one obtains expressions analogous to (23)- 
(25), for the ensemble mean square vorticity components averaged over y. Isotropy again requires 
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that these three mean square vorticity components are equal. The condition that, averaged 
over y. (|uq| 2 ) = (|u>3| 2 ) is satisfied for any function p( A,/x), and the demand that (|a>i| 2 ) = 
(|o> 2 | 2 ) yields the relation 


1 ' 
2 
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(1 -h A) + A/i + A 4- 2 + ~ 



(28) 


Solving (26) and (28) for positive A and p then gives 


A = 


1 

2 


and p 



(29) 


which is enforced by taking 

p(A,^) = 5 (A - |) 5 (^ - -^) (30) 


The remaining function 'Piki, k 3 ) in the covariance (22) will be determined from a consideration 
of the two-point correlation function and the energy spectrum tensor. 

Since the disturbance field is partially homogeneous, (Lumley, 1970) that is, homogeneous in 
the x and z directions and inhomogeneous in the y direction, the two-point correlation function 

is - 

Rij (y, y' , r, t) = Aj d 2 kV(k u k 3 )e tkr ~ ivkH u*(y, k)u,(y',k) (31) 

where r = (x l — x,0 — z), A is a constant proportional to the initial disturbance kinetic 
energy, and V(ki,kz) again plays the role of a probability distribution for the disturbance second- 
moments. In order to obtain explicit expressions for the components of Rij(y , y , r, f), the values 
for the parameters A and /i in (29) must be substituted into the component disturbance mode 
equations (17)— (19)- In the boundary-free case, an exponential time decay with coefficient 
—2 vk? is shown in (4); whereas, in this wall-bounded case, the coefficient —3 is found. This 
is because of the value of A, which is zero in the boundary- free case and 1/2 in this wall-bounded 
case. 

An energy spectrum tensor can also be defined at each fixed y as a Fourier transform in x 
and 2 : 

Eij{y,h,h,t) = AV{kx,kz)u*{y,k)iij{y,k)e~ Zuk7t . (32) 

Once again, it is necessary that the energy spectrum be analytic at k = 0. From (17)— (19), 
this analyticity requires that the distribution 'P(ki,k , i) have a leading-order behavior of k . In 
addition, the finiteness of the initial disturbance energy implies that the spectrum should go 
to zero sufficiently rapidly for large k, specifically at t = 0. These requirements lead to the 
normalized distribution 

HhJ = (33) 

7 r 

where a is related to the variance of the probability distribution. Substituting the probabil- 
ity distribution (33), and the mode disturbance equations, (17)— (19), into the expression for 
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the energy spectrum tensor (32), and integrating, it is found that in the freestreain the mean 
disturbance energy, averaged over one y wavelength, 27 t/\/A k = 2\f2n/k, is 


K{t) = ^ J d 2 k J dyEa{y,ki,kz,t) 

= Kq (l + , ( 34 ) 

where Kq is the initial mean disturbance energy, averaged over y. This then identifies A, 
introduced in the two-point correlation equation (31), as A = 4aK 0 /9. 

In the presence of a boundary, the mean disturbance energy decays according to a power law 
with exponent —2, in contrast to the exponent —5/2 obtained for the boundary-free case. The 
difference between the two cases arises from the fact that for bounded flow only two independent 
wavenumbers are integrated in (34), while the integral in (7), for boundary-free flow, contains 
three independent wavenumbers. 

From Eq. (14) with p = 2, it is seen that for laminar, linear decaying disturbances in the 
presence of a wall, the coefficient in the mean disturbance dissipation rate equation is C e 2 = 1 . 5 . 
In the turbulent case, the experiment of Thomas and Hancock, (1977) with a turbulent Reynolds 
number Re T (= K 2 /ue) of 2000, yielded a decay law near the wall of approximately 1, which 
in turn gave a value for C e 2 of 2. The recent direct numerical simulation of Perot and Moin, 
(1995) with Rer = 137, also displayed a decay law near the wall of approximately 1. 

It may appear surprising that the wall affects the decay rate of linear disturbances infinitely 
far from the wall. The derivation of this result is based on the assumption that the initial dis- 
turbance field is an ensemble consisting of linear superpositions of modes which each satisfy the 
wall boundary condition rq(0) = 0. This boundary condition fixes the phases of the disturbances 
in the direction normal to the wall up to infinity, while the phases in the directions parallel to 
the wall are random. The situation is different for decaying isotropic turbulence, where in the 
initial stages of decay nonlinear interactions among the modes randomize the phases in the y- 
direction outside of a viscous sublayer next to the wall. Even when the turbulence has decayed 
to a stage where the nonlinear terms can be neglected, the phases will remain random far from 
the wall. For this reason, turbulence in the final stage of decay is unaffected by a wall outside 
of the viscous sublayer. 

It is finally necessary to consider the condition for the validity of the linear approximation for 
decaying disturbances. Since there is no mean flow U, one cannot impose the obvious condition 
Jpl 1. In the absence of mean flow, the disturbance velocity obeys the Navier-Stokes equation 


d'Ut dui dp d 2 Ui 

dt JtUj dx 3 dxi +U dxjdxj 


(35) 


The nonlinear terms can 
forces, that is 


be neglected when the inertial forces are 


du t 



«: 


d 2 Ui 

dxjdxj 


much smaller 


than the viscous 
(36) 
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The scale for the amplitude of the disturbance velocity is set by the initial disturbance kinetic 
energy, u t ~ Kq . The length scale L is determined from the probability distribution for distur- 
bance second moments, the standard deviation fixing a wave vector scale L~ l . For a probability 
distribution of the form k 2 e~ ak \ L ~ a^, and the condition (36) gives a Reynolds number 

criterion 

v ^ 0 < 1. (37) 

V 

Although the parameter a(= L 2 ) appearing in the probability distributions (10) and (33) does 
not affect the value of the destruction of dissipation coefficient C e 2 , it does appear in the crucial 
inequality (37) which determines the validity of the linear approximation. This inequality can 
be expressed in terms of the disturbance Reynolds number 

Re d = ^ (38) 

ve 0 

where the initial disturbance dissipation rate £o is proportional, in both the boundary-free and 
wall-bounded case, to vKq/cl. Thus, Red ~ ciKq/ iP " , and the condition for the validity of the 
linear approximation for decaying disturbances can be expressed in the form Red 1. This 
inequality gives a quantitative definition of linear disturbances in the laminar regime. 

4 Disturbance Fields in Homogenous Shear Flow 

The theoretical development here parallels that of Townsend (1970), who described the 
structure of turbulence in a free shear flow as a product of the finite distortion of parcels of 
turbulent fluid. This same RDT framework was used by Hunt and Carruthers (1990) (neglecting 
viscosity) to gain further insight into the structure of turbulence in shear flows and to describe 
some inhomogeneous flows. In this study, the relationship of RDT with linear stability theory 
(Speziale, Abib and Blaisdell, 1996; Salhi, Cambon and Speziale, 1997) is expanded to include 
an analysis of the transport equations for the ensemble-averaged disturbance kinetic energy and 
the disturbance energy dissipation rate. As in the turbulence case, such model equations require 
closure through the specification of closure constants. In the two-equation K — e formulation, 
only the disturbance dissipation rate equation contains modeled terms which have unknown 
closure constants. The homogeneous shear flow is used as a calibration flow for the production- 
of-dissipation and destruction-of-dissipation terms. 

In terms of dimensional coordinates (xi,X 2 ,X 3 ) in a fixed frame, the mean velocity is given 
by U\ = Si 2 , U 2 = U$ = 0, where the mean shear S = constant, and the disturbance velocity 
and pressure are denoted by Uj and p, respectively. For the problem of homogeneous shear, it 
is convenient to work (Rogallo, 1984) in a moving frame (xi,X 2^3)> in which the local mean 
velocity is zero and the coordinates are given by 

x\ = x\ — ( Si )x 2 
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X 2 = %2 

x 3 = x 3 , (39) 


where t is time. The mean shear S and the viscosity u determine the time scale T = S 1 , length 
scale L = V velocity scale U = V vS, and pressure scale P = U 2 = vS for the flow. This 
leads to the introduction of dimensionless variables x } = Xj/L (fixed frame), x' } = Xj/L (moving 
frame), t = t/T, Uj = iij/U, and p = p/P. 

In the moving frame, the disturbances obey the incompressibility condition 


du\ dU 2 dU 2 dU 3 _ 

dx\ 0x[ dx 2 dx' 3 


(40) 


and the linearized Navier-Stokes equations given by 


du i 
~dt 


+ U 2 + 


du 2 

~dt 


t 


dp 


+ 


du 3 

dt 


dp 
dx i 
dp 


dx\ dx 


+ 


dp 

dx'i 


0 +' 2 ) 

(>+‘ 2 ) 

(■+ f2 ) 


d 2 rti _ d 2 u i d 2 ui cPui 

dx' 2 dx\ dx' 2 dx 2 dx 3 

d 2 u 2 _ dPu2 (Pu2 d 2 u 2 

dx' 2 dx\dx' 2 dx 2 dx 3 

<9 2 u 3 _ o. d 2 Uj d 2 u 3 d 2 u 3 

dx' 2 dx\dx' 2 dx' 2 dx 3 


(41) 


(42) 

(43) 


The linear approximation is based on the assumption that the nonlinear terms Ujgfj-, in the 
disturbance momentum equations, can be neglected with respect to both the viscous term uV 2 u t 
and the mean shear term u ] dU l /dxj. This approximation remains valid as long as the amplitude 
of the disturbance velocity remains small compared with the velocity scale \/uS. 

An arbitrary solution of this system in the moving frame can be written as a linear super- 


position 


uj(x',t) = f d 3 k 

p(x', t) = J d 3 k p(k, t)e tk x , (44) 


where the mode amplitudes fj( k, t) and p(k, t) satisfy 


kifi + k' 2 f 2 + k 3 f 3 = 0 (45) 

f\ + f2 + ik\p — —K?f i (46) 

f 2 + ik' 2 p = —K. 2 f 2 (47) 

h + ihp = -£ 2 / 3 , (48) 


with k' 2 = k 2 - tki, K. 2 = k 2 + k 2 + fc 3 , and the wave numbers (Aq, k 2 ,k 3 ) are constant in time. 
The solution for the velocity disturbance modes can be written in the form 

/ t (k, t) = Afy(M)/ j (k,0)e-« (k - t) , (49) 
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where the exponent 


<y(k, t) —j K? c It = -k 2 t 3 — k\tc2 t 2 + {k l k l )t (50) 

arises from the viscous damping of all modes. The transfer matrix M has non-vanishing com- 
ponents 

Mu = ^33 = 1 

k 2 

M 22 = ^2 

„ , kjk 2 (0 - 0 O ) , k\t{k\ - k 2 k 2 + kj) 

Mn = kJ? H + 

fc 3 fc 2 (0 o -0) , kik 3 t(kf - k 2 k 2 + k\) f ^ 

32 = — *£ ’ (51) 

with i / u \ 

k H = (ki + kt)\ 00 = tan- 1 (g-) , 0 = tan- 1 (g-J (52) 

The initial mode amplitudes /,(k,0) are elements of an ensemble with mean zero and a co- 
variance (/j(k, 0)/j(l, 0)). Assuming that the initial disturbance field is homogeneous, isotropic, 
and has finite energy Kq, the initial covariance is given by the same expression as above for the 
case of zero mean shear with no boundary 

</*(M)/;( 1,0)) - 5 3 (k + 1) fa - ^) KoVQt), (53) 

where 2 

p (k> = i (£) ^- ak ’ ( 54 > 

is the probability distribution for disturbance second-moments, and the parameter a, related to 
the variance of the distribution, will be given a physical interpretation in what follows. This 
probability distribution is consistent with a k 4 low wave number behavior for the energy spectrum 
E(k) (Batchelor and Proudman, 1956) and, in the case of linear disturbances in zero-mean-shear 
flow, leads to a decay law for the mean disturbance kinetic energy. 

The disturbance mode covariance is related to the energy spectrum tensor by 

(fi(k,t)fj(l,t)) = $y(M)* 3 (k + l), (55) 

so it follows from Eq.(49) that 

Eij( k, t) = Mu( k, t)M jm ( k, i)^oim(k)e- 2 ^ k - f ), (56) 

where, as a result of (53), 

Eoijik) = Eij{ k, 0) = (Sij - K 0 V(k) (57) 


13 



is the initial energy spectrum tensor. 

The energy spectrum tensor is the fundamental quantity from which all disturbance second- 
moments are derived by integration over wave number space. Of particular interest are the 
disturbance stress tensor 

Tij(t) = I d 3 k£ tJ (M) , (58) 

the disturbance kinetic energy K(t) — Tji/2, and the disturbance dissipation rate 

e{t) = J d 3 k K? E lt {k,t). (59) 

A physical interpretation can be given to the parameter a, which appears in the probability 
distribution (54), by performing the integral in (59) at time zero to obtain eo — 5Ko/a. In terms 
of the dimensional disturbance kinetic energy K = vSK, and the dimensional dissipation rate 
e = vS 2 e the parameter a is given by 

a = 5 = 5t7o (60) 

It will be shown that this quantity, a ratio of the disturbance time scale to the mean flow time 
scale, regulates the subsequent growth or decay of linear disturbances in homogeneous shear 
flow. In the turbulence case, the variable rj = SK /i also plays a critical role. Jongen and Gatski 
(1998) have shown that for turbulent homogeneous shear flow, t) reaches an equilibrium value 
which can be analytically connected to the variation of the ratio of kinetic energy production to 
dissipation rate. 


5 Modeled Disturbance Transport Equations 


The ensemble-mean kinetic energy and dissipation rate for linear disturbances satisfy evolu- 
tion equations which can be derived from the linearized Navier-Stokes equations (in the fixed 
frame). In this homogeneous flow, gradients of all ensemble-mean quantities vanish, and the 
(dimensionless) kinetic energy and dissipation rate equations reduce to 

K = — t 12 -e = V K -e, (61) 

£ = V £ - V e , (62) 


where 7 7 a'(= -7*12) is the kinetic energy production, 

v 0 (/ ^l d*h \ / dii 1 du 2 \ \ 

He ~ \\dx l dx 2 r\dx l dx l / ) 

= -2 J d 3 k \k x k' 2 E u + K 2 E n ] 

is the production-of-dissipation, and 


V £ =2 


d 2 Ui d 2 
dxjdxk dx 


%) = 2 J’ 


d 3 k/C 4 Eu 


(63) 


(64) 
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is the destruction-of-dissipation. From Eqs. (61) and (62), it is seen that only the disturbance 
dissipation rate equation requires closure models. In the standard K — £ formulation, the 
production-of-dissipation is modeled by 



Ve = -C e i— Ti2, 

(65) 

which gives 

C] - v K 

£l e £T l2 £ Vk' 

(66) 

for the production-of-dissipation 

coefficient. The destruction-of-dissipation is modeled by 



£ 

II 

to 

(67) 

which yields 

_ K _ 

C £ 2 = ~V> £ 

£ l 

(68) 


for the destruction-of-dissipation coefficient. 

Recall that the components E l} of the energy spectrum tensor are completely defined by 
Eqs. (49) through (57). With the solution of these equations, all the quantities of interest, 
including V e and can be obtained by carrying out the integration in Eqs. (58), (59), (63), 
and (64). These three-dimensional integrations were initially performed numerically using cubic 
and spherical grids. It was found that although the convergence was good at small times, at 
times greater than about 105' 1 , the results became strongly dependent on the grid resolution 
used (even for grids with as many as 250 million points). The numerical evaluation of the 
integrals in Eqs. (58), (59), (63), and (64) is therefore greatly simplified, and the accuracy 
significantly enhanced, by performing the radial integration analytically. 

It is convenient to introduce spherical coordinates (A;,0, </>), in wave- vector space, such that 

k\ = k sin 9 cos 0, k 2 = k cos 6 , £3 = k sin 9 sin <fi, 

and ktf — fcsin/9. For the chosen probability distribution (54) the radial integrals in Eqs. (58), 
(59), (63), and (64) are Gaussian and can easily be carried out analytically. The ensemble- 
averaged disturbance stress components are then given by 

is t"K r2n 5 

nAt) — / d9sin9 d4f-*eij, (69) 

47t J o Jo 

where the angular functions / = ), and e tj = eij{6,4>,t) are defined in the Appendix in 

Eq. (Al) and Eqs. (A2) to (A6), respectively. This yields for the ensemble-averaged kinetic 
energy 

K r n r 2 * 5 

K{t) — — ^ / d9sin9 d<t>f-2e lU (70) 

o7T Jo Jo 

and for the ensemble-averaged disturbance dissipation rate 

e{t) = f dO s'm 9 f d(j)f~^eup (71) 

87T7?o Jo Jo 
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where p = K, 2 jk 2 so that 

p = p(Q , fa t) = 1 - t sin 20 cos <f> + t 2 sin 2 0 cos 2 <f> (72) 

The resulting expressions for the production-of-dissipation and destruction-of-dissipation 

K C n /•27 r 7 

= — / dQs'mO d<f>f~*[e ti s + enp], (73) 

47r7/ 0 y 0 y 0 

£> e (£) = f dO sin 0 f d(j)f~^e zl p 2 , (74) 

407t?7q Jo Jo 

respectively, where 

s = s(0, <f>,t) = - sin 20 cos <j> — £sin 2 0cos 2 <f> (75) 

Convergence tests show that these integrals had to be evaluated on grids with spacing of 1/16 
of a degree. With the evaluation of the integrals in Eqs. (69) - (71), Eqs. (73) and (74), it is 
now possible to analyze the ensemble-averaged transport equations for the disturbance kinetic 
energy and dissipation rate. Of particular interest is the dissipation rate equation Eq. (62) 
which requires closure. It will be shown in the next section that closure models for T* e and 
can be determined through an analysis of the long-time behavior of the disturbance correlations. 

6 Results 

From an examination of Eqs. (69) - (74), it is apparent that the only factor containing the 
parameter Tjo is /(0, given in Eq. (Al). This factor, which appears in all the ensemble- 
averaged correlations and which arises from viscous effects, governs the large-time behavior of 
all correlations causing them to ultimately decay. At large times and generic angles sin 0 cos (j) ~ 
0(1), the i 3 -term in / dominates, leading to negligibly small values of the integrand. However, 
for sin 0 cos 0 = k x /k small, / is proportional to t at large times. The t 2 - and t 3 -terms in f will 
also be 0(t) when (0,^) lies in a band of solid angle of width 0(t *) about the k\ = 0 plane. 
At large times, this band of solid angle gives the dominant contribution to the integrals in Eqs. 

(69) -(74). 

Now, consider the disturbance kinetic energy (70). While / 2 decays at all times, the factor 

e ti (see Eq. (A5)) grows, indicating that the evolution of the disturbance energy is governed by 
a relative balance of the two factors. Figure 1 shows the evolution of the disturbance kinetic 
energy for different initial values of rj 0 . For a small value of rjo (= 0-6), the factor f ~ 5/ 2 dominates 
and the energy decays monotonically in time. A transition between this monotonic decay and 
substantial growth of the disturbance energy occurs for values of r/ 0 between 1.2, where the 
energy reaches a minimum and grows briefly, and 2.4, where the energy grows to a magnitude 
just equal to its initial value before decreasing. For larger values of rjo (> 3.0), the viscous decay 
factor is further suppressed and, after a brief initial period of decay, the energy then grows for a 
period of time before finally decaying. These larger values of tjo can be compared to results from 
the turbulent case. The initial value of r/o = 6 is a moderate shear case which corresponds to 
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the equilibrium value for 77 in turbulent homogeneous shear flow. The temporal evolution shown 
in Fig. 1 for this value of 770 contrasts with the exponential energetic growth (Tavoularis, 1985) 
characteristic of the turbulent case. Here the disturbance kinetic energy grows to a little beyond 
three times its initial value. Therefore, if the initial disturbance kinetic energy is considerably 
less than the natural energy scale vS, the linear approximation should remain valid for all times. 
Calculations have also been carried out for 77 0 = 12 and 18. The evolution of K with time for 
these other 770 values (not shown) is very similar to that for r/o = 6 . The only difference is that, 
as expected, the maximum value of K reached before the final decay is larger than for 770 = 6 , 
for 770 = 12, Kmax = 6.95Ko, and for 770 = 18, K max = 10.82 K 0 . The initial value of t ? 0 = 18 
is a high-shear case which is very close to the value used in the DNS study of Lee et al. (1990) 
(770 ss 17). Using the numerical simulation results, Lee et al. also compared with RDT and found 
very good agreement over the time interval examined (0 < t < 12). (The quantitative basis of 
comparison between the DNS and RDT were the Reynolds stress anisotropy components which 
are discussed below and shown in Fig. 4. The high-shear limiting value 770 -7 00 corresponds to 
the case traditionally considered in RDT calculations (Hunt and Carruthers, 1990; Rogers, 1991), 
where viscosity is neglected. The absence of viscosity allows the energy to grow monotomcally, 
rendering at some finite time the linear approximation invalid. 

Figure 2 shows qualitatively similar trends for the disturbance dissipation rate e{t) at the 
corresponding values of 770 = 0.6 and 6.0. For the small initial value 770 0.6, the dissipation 

rate also decays monotonically. At the moderate shear rate value 770 = 6, e(f) grows after a brief 
initial period of decay, reaches a peak and then decays at a more rapid rate than the disturbance 
kinetic energy. The more rapid decay of e(t) compared with K{t), at large times, comes about 
because appears in the integral (71) for the dissipation rate, while the lower power / 2 
appears in (70) for the kinetic energy. Similar behavior was found for the higher shear cases of 
770 = 12 and 18. 

Equation (61) describes the evolution of the disturbance kinetic energy with time. The 
equation shows that for Vk! e ratios less than one, the kinetic energy must be a decaying function 
of time. Figure 3 shows that for the very weak shear initial condition of 770 — 0.6, the VkI £ ratio 
never reaches unity confirming the continuous energetic decay of the disturbances. However, for 
the initial condition 770 = 6 , 12 and 18, the V K /e ratio exceeds unity for a period of time and 
then decays below unity. This latter behavior is reflected in the disturbance energy evolution 
shown in Fig. 1 for 770 = 6 where the disturbance energy grows up to a time of approximately 25 
(Ptc/e > 1) and then decays (P/c/e < 1) at larger times. This same qualitative behavior holds 
for the higher shear rate cases of 770 of 12 and 18 as well. 

Since the main focus here is on disturbances which can eventually grow (and lead to tur- 
bulence), the results for the initial value of 770 = 0.6 are not of interest. In addition, while 
this weak shear case yields decaying disturbances, the validity of neglecting the nonlinear terms 
does come into question. As Cambon and Scott (1999) have pointed out, the contribution of 
the nonlinear term, while small, may have a cumulative effect on the dynamics over a period 
of time. In addition, the linearized, weak shear case may be inconsistent because the product 
term containing the weak mean shear is retained while the nonlinear terms, possibly of the same 
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magnitude, are omitted. These considerations dictate that subsequent results will focus on the 
moderate- and high-shear cases. 

Figure 4 shows the evolution of the ensemble-mean stress anisotropies 




2 K(t) 



(76) 


for the initial condition 770 = 6.0. It is seen that initially isotropic linear disturbances are 
distorted by mean shear, rapidly becoming anisotropic. The stress anisotropies asymptote to 
the values b n = 0.631, b 22 = -0.333, 633 = -0.298, and b l2 = -0.008, which is consistent with 
the turbulent DNS (and RDT) results (Lee, et a/, 1990) at high-shear rates. In the current 
study, the bn, b 22 , and 633 anisotropy components are computed individually. The resulting 
trace (bn = 0) is then used as a measure of the accuracy of the simulation which, for the results 
presented here, was satisfied to C?(10 -15 ). 

Recall that initially, an isotropic distribution of the disturbance energy is assumed. Thus, 
the structure of the disturbance field is significantly altered by the imposition of the mean 
shear. This alteration is easily seen in the mapping of the anisotropy invariants //(= —bijbij/2) 
and ///(= b lk b kl bij/3) shown in Fig. 5. The figure shows that after t « 2, the (realizable) 
disturbance field migrates toward the two-component (2C) state. After t ss 10 the disturbance 
field is 2C and monotonically evolves towards a one-component (1C) state. Also shown in Fig. 
5 is the invariant map trajectory of the DNS results (Lee, et al, 1990). As can be seen, at early 
times ( t < 2) both the DNS trajectory and that from the present calculations, with r/ 0 = 6 and 
18, are in phase, and at later times, even though the DNS evolution lags behind the the current 
results, all of these evolve toward the same 1C state. At large times, the system has evolved 
to what has been termed a “statistical eigensolution” of the moments (Hunt and Carruthers, 
1990). That is, if the disturbances were to have these forms initially, the disturbance field would 
change little under a distortion. While the results shown in Fig. 5 are for the 77 0 = 6 and 18 
cases, the same trend toward a 1C limit was found to hold for the case of 770 — 0.6. Thus, for 
linear disturbances the existence of the statistical eigensolution is not limited to the high-shear 
case. These disturbance field results are in sharp contrast to the equilibrium state reached for 
the turbulent case where the fixed point invariant values are IIIqo ~ 0.0043 and II<x> ~ 0.064. 
These values are far removed from the 2C and 1C states exhibited by the linear disturbance 


field. 

The results up to this point have ascribed a structure to the linear disturbance field in 
homogeneous shear which differs from the turbulent homogeneous shear case. In addition, the 
time scale ratio 77 does not evolve (see Fig. 6) to a constant value, as in the turbulent case, but 
displays linear algebraic growth at large times. For a particular homogeneous flow ( Jongen and 
Gatski, 1998) , the dynamical behavior of the system can be described by the time scale ratio 
77, and the production-to-dissipation rate ratio ’P/c/ £ discussed previously. 

The consequences of these results for developing a closure model for the disturbance dissi- 
pation rate equation (62) can be seen from the evolution equation for 77 given by 

f) = ~(C £l -1)^ + ((7*2-1), (77) 
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As seen from Fig. 6, 77 displays a linear behavior at large times so that 77 (= 0.34, for 770 = 6) is 
constant. The large time value of 77 is nonzero and depends weakly on t/o- This is in contrast to 
the turbulent case where at large times 77 = 0. Since the production-to-dissipation rate ratio is 
not constant (see Fig. 3), Eq. (77) suggests that the closure coefficients C £i and C e 2 may also be 
functions of time. Figure 7 shows the evolution of the disturbance dissipation rate coefficients. 
Results for three initial values of 770 are shown which suggest that a limiting range of values at 
large times (and large 770) can be reached. For C £ 1 there is a modest change in the value at 
t = 50 with changes in 770: for 770 = 6, C £ \ = 2.23, for 770 = 12, C £ 1 = 2.08, and for 770 = 18, 
Cel = 2.01. For C e 2, the values range about ±3% from a mean of 2.46. Although it was not 
possible to determine analytically the limiting values for C £ 1 and C £ 2, numerical simulations for 
large values of tjq at large times showed that the limiting values were C e \ ~ 2.0 and C e 2 ~ 2.5. 
These values are in contrast to the usual values obtained for turbulent closure models where 
C £ 1 ~ 1.45 and C £ 2 lies in the range of 1.83 to 1.92. 

In the turbulent case, the destruction-of-dissipation coefficient was deduced from an analysis 
of the decay of isotropic turbulence (e.g. Comte-Bellot and Corrsin, 1971). In the analysis of the 
decay of homogeneous, isotropic linear disturbances, above, a value for C £ 2 of 1.4 was determined 
in a boundary-free case. This is in sharp contrast to the value of 2.5 found in this homogeneous 
shear flow. An explanation for this difference may lie in the structure of the disturbance field. 
In the previous decay studies, the disturbance field was isotropic, and remained isotropic in the 
absence of mean shear. However, as the results shown here indicate, the imposition of mean 
shear quickly produces an anisotropic field (see Fig. 4) which eventually drives the flow close to 
a 1C state (Fig. 5). Thus, a more relevant flow with which to form a comparison may be the 
decay of anisotropic turbulence. Dakos and Gibson (1990) studied such a flow and deduced from 
the decay of the turbulence the destruction-of-dissipation rate coefficient C £ 2 . Their anisotropic 
decay data yielded a value for C £ 2 of 2.18 which is relatively close to the value of 2.5 found here 
for C £ 2 - Of course, in spite of fundamental differences in the dynamics associated with each 
study, it is interesting to note that the introduction of anisotropy in either disturbance field 
(linear or turbulent) significantly increases the value of the coefficient C £ 2 . 

The value of the production-of-dissipation rate coefficient C £ i used in turbulence modeling 
is determined from a consistency criterion found from homogeneous shear flow at equilibrium. 
As alluded to earlier, the time scale ratio 77 is constant (77 = 0) in the turbulence case so that 
Eq. (77) yields the well-known relation 


(Vk\ = < 7 . 2-1 
V £ ) 00 <7,i- r 


(78) 


With the value of 7>k/£ fixed from experiments (for example), and C £2 fixed from the decay of 
isotropic turbulence, the value of C £ 1 can be determined. Unfortunately, this is not the case 
here and, as Fig. 7 shows, the value of C £ \ varies with time. As noted previously, it was not 
possible to analytically determine the limiting values for C £ 1; however, simulations were run at 
large values of 770 and for long times to try to determine its asymptotic limit. These simulations 
showed that C £ \ approached 2.0 at long times which is in contrast to the value of C £ \ ~ 1.45 for 
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the turbulent case. Thus, for both the production-of and destruction-of-dissipation rate models 
the closure coefficients are significantly larger than the corresponding turbulent values. 

7 Summary and Conclusions 

This study has introduced an approach to transition modeling in which deterministic solu- 
tions of the linearized Navier-Stokes equation are combined with a probability distribution that 
accounts for the uncertainty in initial conditions to obtain mean transport equations for distur- 
bances in the laminar regime. This approach was applied to the study of decaying disturbances 
in zero mean-shear flows. The analysis has shown that the linear disturbance field in laminar 
flow decays at rates which differ from that of the turbulence field. In the boundary-free case, the 
disturbance kinetic energy decays at a much faster rate of 5/2 than the corresponding turbulent 
kinetic energy decay rate of ~ 1.25. This yields a destruction-of-dissipation rate coefficient C e 2 
in the mean dissipation rate transport equation of 1.4 for the linear disturbances and ~ 1.8 for 
the turbulence. In the case with a wall boundary, the disturbance kinetic energy decayed at a 
rate of 2, which is faster than the corresponding turbulent kinetic energy rate of 1.0. This yields 
a destruction-of-dissipation rate coefficient C e2 of 1.5 for the linear disturbances and « 2 for the 
turbulence. These results are summarized in Table 1 which clearly shows that the wall has the 
effect of reducing the decay rate of the disturbance energy in both the laminar and turbulent 
regimes. 

Closure models for the production-of-dissipation rate and destruction-of-dissipation rate in 
a linearized disturbance transport equation for the dissipation rate have been developed. In 
this study of mean homogeneous shear flow deterministic solutions of the linearized Navier- 
Stokes equation are combined with a probability distribution that accounts for the uncertainty 
in initial conditions, to obtain mean transport equations for disturbance correlations in the 
laminar regime. 

The temporal evolution of the disturbance kinetic energy and the disturbance dissipation 
rate were shown to depend on the magnitude of the initial value of the time scale ratio r/. The 
temporal evolution for both the kinetic energy and dissipation rate was monotonically decaying 
for values of 770 < 1.2. For larger values of 77 0 , both quantities initially decayed, then grew 
with time, and then subsequently decayed at large times. This is in contrast to the turbulent, 
homogeneous shear case where the quantities displayed exponential growth. 

Even though the temporal behavior of the disturbance kinetic energy differed from the tur- 
bulence case, the disturbance stress anisotropies did reach a steady-state which corresponded 
to a “statistical eigensolution” of the problem. An analysis of the anisotropy invariant map 
shows that this eigensolution corresponded to the 1C limit for the disturbance field, which is in 
contrast to the structure shown by turbulence at equilibrium (b^ = 0). 

In the turbulent, homogeneous shear case, the time scale ratio 77 is a fixed point of the 
system; whereas, in the linear disturbance case it was found to have an 0(t) algebraic growth. 
As in the turbulent case, the evolution equation for 77 yields a relation between the Vk/z ratio 
and the closure coefficients in the disturbance dissipation rate equation. Unlike the turbulence 
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case, however, the r] and Vk/z variables are time dependent which leads to time dependent 
behavior for dissipation rate coefficients C £ \ and C £ 2. Nevertheless, limiting values (for large 770 
and at long times) were obtained which yielded values for C e \ and C € 2 that were larger than the 
corresponding turbulent values. 
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Appendix 


The angular functions /( 0 ,<M), and e tj = e y ( 0 , 0 , t) appear in Eqs. ( 69 ) to ( 74 ), 

and are 

given by 

2 q(k,t) 

/(».*«) -‘♦aaV 

= 1 + — ft - -t 2 sin20cos</> + ^< 3 sin 2 0cos 2 <jA , 
brio \ 2 a / 

(Al) 

and 



en(0,^,t) = (cos 2 0 + sin 2 0 sin 2 4 >) 



+ (2 cot 0 sin 2 </>) a + (esc 2 0 sin 2 (f> tan 2 <^j a 2 
- (sin 20 cos 3 4 > + 2 a sin 2 </> cos <pj ^ 



/ r ^\ 2 

+ sin 2 0cos 4 (/> (— J , 

(A2) 


sin 2 0 

e 22(0> — ^2 ’ 

(A 3 ) 


e 33 (0, <£, t) = (cos 2 0 + sin 2 0 cos 2 </>) 



— (2 cot 0 sin 2 a + (esc 2 0 sin 2 a 2 



vt 

— sin 2 (f> cos <p (sin 20 — 2a) — 

/ r A 2 

+ sin 2 0 cos 2 0 sin 2 ? 

(A 4 ) 


eu( 0 ,<f>,t) = 1 + cos 2 0 + 


sin 2 0 N 


+ 


(esc 2 0 tan 2 <f>^ a 2 


rt n „ 9 , / rt \ 2 

- (sin 20 cos 0 ) — + sin 9 cos <f> J , 


(sin 20 cos cf> + 2a sin (f) tan 4 >) 
ei2{0,(p,t) = ^ 

/ sin 2 0 cos 2 <f>\ rt 
+ \ P j p' 

where p = p( 9 , 4 >, t) is given in Eq. ( 72 ), and 

a( 9 ,<fr,t) = tan -1 (cot0) — tan -1 (cot0 — tcosrfi) 
r( 9 ,<f>,t ) = - cos 20 + ^ sin20cos</>. 


(A 5 ) 


(A6) 

(A 7 ) 

(A8) 
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Case 

p 

c £2 

Laminar-No Wall 

2.5 

1.4 

Laminar- Wall 

2.0 

1.5 

Turbulent-No Wall 

» 1.10 

Si 1.9 

Turbulent-Wall 

» 1.0 

ss 2.0 


Table 1 . Decay exponent p and model coefficients C £ 2 for laminar and turbulent cases. 
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Figure 5: Invariant mapping of disturbance anisotropy tensor at various times for initial con- 
ditions: r] o = 6 (times as labeled); r) o = 18; A, rjo « 17 from DNS (Lee, et al, 1990) at 
t = 2,4,6,8,10,12. 





